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Abstract 

Generalized Bose-Einstein symmetry requires that J + I = even for two-pion angular states of spin J 
and total isospin /. We show that the symmetry predicts three linearly independent constraints on partial 
wave intensities with even spin for ir^p — > vr^vr+n, ir^p vr^vr^n and vr+p — > Tr+vr+n. Available data 
violate all three constraints for S-, D^-, D^- and Z?^-partial waves. The violations of the symmetry imply 
a presence of the symmetry violating contributions in transversity amplitudes in Tr~p — > '!T~'!T^n and predict 
quantum entanglement of 7r~7r+ isospin states which is excluded by the symmetry. We derive approximate 
lower and upper bounds on entanglement amplitudes \as\ and \aA\- The bounds provide a clear evidence for 
entanglement of 7r~7r+ isospin states below 840 MeV and suggest the entanglement at higher dipion masses. 
The small values of \as\ ^ 0.10 — 0.20 below 840 MeV explain the puzzling differences between the S'-wave 
intensities in Tr~p — > 7r~7r+n and n^p —^ Tr^ir^n and reveal a suppression of isospin / = 0, 2 contribution in 
the S'-wave amplitudes in Tr~p — > vr^vr+n. The large isospin I — 1 contribution of p'^{770) to both S- and 
P-wave amplitudes is due to large entanglement amplitude |a^| ~ 0.98 — 0.99. These findings confirm the 
predictions of a model of non- unitary dynamics of the pion creation processes arising from a CPT violating 
interaction of these processes with a quantum environment. 
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I. INTRODUCTION. 



Measurements of pion creation processes 7r~p — > 7r~7r+n [H, S, S, 0, [H] and vr+n — > 7r+7r~p 0] on 
polarized targets enable model independent determination of production amplitudes. Evidence for 
a narrow rho-like resonance in the S-wave transversity amplitudes was found in amplitude analyses 
of CERN data on Tr~p — > 7r~7r^n at 17.2 GeV/c at small momentum transfers 0,(3)13)0) l3) 11: l3] 
as well as at large t \4^, in analysis of ITEP data on the same reaction at 1.78 GeV/c at small 
t and in the analyses of CERN data on vr+n —>■ 7r"^7r^p at 5.98 and 11.85 GeV/c at larger 
momentum transfers 0, [l^]- The S- and P-wave subsystem of the reduced density matrix 
measured on transversely polarized target is analytically solvable at any dipion mass The 
moduli of the S'-wave transversity amplitudes \Sr\'^ are given by a simple expression 



\St\ = ai,T + a2,T - 3|Lr| (1-1) 

where t = u,d are target nucleon transverse spins "up" and "down" relative to the scattering 
plane. The terms ai^j- + 02, r are expressed in terms of measured spin density matrix elements and 
exhibit a clear p^{770) resonant peak [l^ that survives the subtraction of the p^{770) peak in the 
P-wave amplitude |i^Tp- The rho-like resonance present in the S'-wave amplitudes - previously 
referred to as (t(750) resonance - must be interpreted as p^{770) resonance. We conclude that the 
data on polarized targets reveal p^{770) — /o(980) mixing in the S'-wave amplitudes which appears 
to violate Lorentz symmetry and Generalized Bose-Einstein symmetry for pion isospin multiplet. 

The origin of p^{770) — /o(980) mixing is related to another facet of the CERN data on polarized 
target. In the previous work [[15|] we showed that the data at large momentum transfers t provide 
evidence for evolution of pure initial states to mixed final states in 7r~p — > 7r~7r+n. In quantum 
theory such non-unitary evolution occurs in open quantum systems S interacting with a quantum 
environment E. While the system S and the environment E undergo a unitary co-evolution 
Pf{S,E) = Upi{S, E)U~^ , the evolution of the reduced density matrix Pf{S) = TrEpAS. E) = 
£{pi{S)) of the system S is non-unitary and is given by Kraus representation [l^, 17. Ha. Eil 



Pf{S) =£{pi{S)) = ^^PmnS^„/)i(S)S+^ (1.2) 

l m,n 

where Sim =< G£\U\em > and |e^ > are interacting degrees of freedom of the environment. The 
initial state of the environment has a general form 

Pi{E) = '^Pmn\em >< Enl (1-3) 

We showed in [3] that Kraus representation leaves invariant the formalism used in data analyses 
provided the co-evolution with the environment conserves P-parity and quantum numbers of the 
environment. On general grounds 

0) 

EtI there are four interacting degrees of freedom of the 
environment in n^p vr^Tr+n. The Kraus representation then has a diagonal form 

4 4 

Pf(^) = ^PuSuPi{S)Sjf^ = ^pepfii) (1.4) 

The question now arises what are the quantum states of the environment in tt~p — > 7r~7r^n. 
In [l^ we associate with the two solutions for transversity amplitudes Au{i) and Aii{j) two 
one-qubit states \i > and |j > where i,j = 1,2 and identify the quantum states \e£ > with four 
two-qubit states \i > \j >. The final state Pf{S) is mixed state of four solutions Pf{ij) with 
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probabilities pi = pij and ^ pij = 1. The measured final state Pf{iT TT~^n) is mixed state even 

when the initial states Pi{TT~p) are pure states. According to Wald Theorem [20'], the interaction 
giving rise to such a non-unitary evolution from pure states to mixed states must violate CPT 
symmetry. The violation of the CPT symmetry implies that such interaction violates also Lorentz 



symmetry 21[ . 



In our previous work [l^ we have formulated a model of non-unitary dynamics of 
p°(770) — /o(980) mixing. According to the model, the pion creation process iT~p ■K~'K^n pro- 
ceeds in three stages. In the first stage resonant and non-resonant (isospin 2) qq (qqqq) modes with 
spin K, helicity /i and isospin Ik are formed in a process ir'p — > qqn that preserves Lorentz and 
CPT symmetry of QCD dynamics. In the second stage the resonant modes propagate with Breit- 
Wigner amplitudes. As the modes propagate they interact in the last stage with the environment 
states \i > \ j >. This CPT violating interaction has two effects. First, the transversity amplitudes 
Au and split into two solutions each Au{i),i = 1,2 and Ad{j),j = 1,2 with the solution labels 
now re-interpreted as quantum numbers of the environment. Second, the interaction induces 
isospin conserving transitions between qq states of spin K, helicity fi and isospin Ix into two-pion 
states with spin J, helicity A and isospin /(7r-7r+) = Ik- As a result, /9°(770) appears in the 



5- wave amplitudes and /o(980) appears in the P-wave amplitudes, as observed experimentally jl3l |. 



Generalized Bose-Einstein symmetry extends the Bose-Einstein statistics for identical bosons 
to the entire isospin meson multiplet [23]. As the result, two-meson states must have a definite 
symmetry under permutation of the members of the multiplet. In the case of two-pion states 
the Generalized Bose-Einstein symmetry requires that J + I = even where J and / are spin and 
isospin of the two-pion state [23| . 

The environment induced transitions between qq modes and two-pion states violate Generalized 
Bose-Einstein symmetry. Transversity amplitudes A"^^ are then a superposition of Generalized 

Bose-Einstein symmetry conserving and violating amplitudes C^'^ and V^'^ 

A'cii) = «JA,,c^:(o+^JA,,<;:(^) (1.5) 
Ai:}u) = «^A,,c5(j)+^jA,,<';(j) 

where tj is i-channel naturality of the amplitude |15 ]. The self-consistency of (1.5) with angular 
expansion of the total amplitude requires [l3l | 

J=even: ajx,rj = as, i^jx,r, = o-a (1-6) 
J=odd: ajx,n = a^, wja,^ = as 

where la^p + |a^p = 1. The novel feature of the model is its prediction that the produced two-pion 
states are not separable but entangled isospin states 

|^(7r"7r+) >= as\S > +aA\A >= a|7r~7r+ > +6|7r+7r" > (1.7) 

where \S > and \A > are symmetric and antisymmetric 7r~7r^ isospin states, respectively, and 

a = -^{as + aA) (1.8) 



b = j^ias 



a A) 



Only for as = aA = A= are the two-pion states vr tt~^ separable. 
v2 
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The assumption of Generalized Bose-Einstein symmetry implies a relation between transversity 
amplitudes for even spin J in tt~p — > 7r~7r+n, Tr~p — > TT^n^n and vr+p — > vr+vr+n 

= 4';(00) - Vf (1-9) 

The symmetry thus imposes constraints on the observables in these three reactions. The purpose of 
this work is to use the available data to test the validity of the Generalized Bose-Einstein symme- 
try. In addition, we examine the data for evidence of quantum entanglement of 7r~7r"'" isospin states. 

The paper is organized as follows. In Section II. we introduce Generalized Bose-Einstein sym- 
metry of two-pion states in reactions vrA^ ttttN. In Section III. we derive three constraints on 
partial wave intensities with even spin in n^p 7r~7r''~n, 7r~p ir^Tr^n and 7r~^p — > 7r"^7r"^n which 
are the consequence of Generalized Bose-Einstein symmetry. In Section IV. we gather data required 
for the tests of the constraints and present the results of the tests. The available data violate the 
symmetry in all partial waves tested. In Section V. we derive approximate lower and upper bounds 
on entanglement amplitudes \as\ and |oa|. The bounds provide a clear evidence for quantum en- 
tanglement of 7r^7r+ isospin states below 840 MeV and suggest the entanglement at higher dipion 
masses. The small values of |a5| ~ 0.10 — 0.20 below 840 MeV explain the puzzling difference 
between the 5- wave intensities in iT~p — > 7r~'ir~^n and tt~p vr^vr^n and reveal a suppression of 
isospin / = 0, 2 contribution in the 5-wave amplitudes in tt~p — > TT~TT~^n. The large isospin 1 = 1 
contribution of p^{770) to both S- and P-wave amplitudes is due to large entanglement amplitude 
\aA\ ^ 0.98 — 0.99. The paper closes with a summary and concluding remarks in Section VI. 



II. GENERALIZED BOSE-EINSTEIN SYMMETRY OF TWO-PION STATES IN 
ttN ttttN REACTIONS. 



A state consisting of two identical particles must be described by a state vector with definite 
symmetry properties with respect to interchange of the particles. The two-particle state must be 
symmetric if the particles are bosons and antisymmetric if they are fermions. Consider two-body 
angular momentum state \JM; //i/i2 > where J is the total spin and /ii and ^2 are helicities of the 
particles. The properly symmetrized or antisymmetrized physical state is given by [i^ 

Ni\JM- ^i/i2 > +(-l)-^| JM; /ism >) (2.1) 

where the normalization coefficient N = if ^ fi2 and N = ^ if fii = ^2- Bose-Einstein 
statistics requires that two-pion states vr'^vr'^ and 7r+7r+ must have even J. 



Generalized Bose-Einstein symmetry is an extension of Bose-Einstein statistics to all particles 
belonging to an isospin multiplet of isospin Ii which are regarded as 2/i + l charge states of the same 
particle. By incorporating the isospin quantum numbers in the state vector the symmetrization 
properties are extended to the interchange of particles belonging to the same multiplet. Then 
the properly symmetrized or antisymmetrized angular momentum state of total isospin / of two 
particles belonging to the same multiplet Ii is given by [i^ 

\JM-fiifi2]Im > +{-iy+^-'^^'\JM; fi2f^i; Im > (2.2) 

If Hi = (J.2 then J + I — 2Ii must be even. Generalized Bose-Einstein symmetry requires that for 
two-pion state J + / is even. 



5 



Consider separable two-pion state [vr > Ivr"*" >. It can be written in the form 

|7r~7r+ >= > > (2.3) 

where \S > and \A > are symmetric and antisymmetric maximally entangled tt~tt^ isospin states 

|5 > = -^{\tt- > |7r+ > +|^+ > Ivr" >) = — ^(V2|0,0 > +|2,0 >) (2.4) 

\A> = ^(Ivr" > |7r+ > -|7r+ > \tt~ >) = |1,0 > 
V2 

where we used the convention |7r^ >= — 11, +1 > [24]. Transversity amplitudes are matrix elements 

A'^^ =< JAr/l < ^-^+1 < r„|T|Or >= ^Ci^l + (2.5) 

where r and r„ are target and recoil nucleon transversities and stands for incident pion he- 
licity. The transversity amplitudes are equal to Generalized Bose-Einstein symmetry conserving 
amplitudes C'{^^ 

J=even: C/;^ = < J\r]\ < S\ < r„|r|OT > (2.6) 
J=odd: C/;;? = < JAr?| < A\ < r„|r|Or > 

since the Generalized Bose-Einstein symmetry violating amplitudes V^l vanish 

J=even: V^'Jl = < JXrj\ < A\ < T„|r|Or >= (2.7) 
J=odd: V('^ = < JXr]\ < S\ < r„|r|Or >= 

The two-pion states Ivr^Tr*^ > and |7r"'~7r^ > satisfy the constraint J + I = even since they correspond 
to even / 



vr vr > 



-^(|0,0>-V2|2,0>) (2.^ 



|7r+7r+ > = |2,2 > 

Generalized Bose-Einstein symmetry is assumed to be a symmetry of strong interactions. It is 
violated by electromagnetic interactions. The observation of /?°(770) — /o(980) mixing suggests it 
is violated also in the interactions of the pion creation process with the quantum environment. 



III. CONSTRAINTS ON PARTIAL WAVE INTENSITIES FROM GENERALIZED 
BOSE-EINSTEIN SYMMETRY. 

We first consider S'-wave amplitudes Sr{ — h), S't(OO) and S't-(++) in three measured processes 
7r~p — > 7r~7r^n, 7r~p n^TT^n and vr+p vr+vr+n 



Sr{- + ) 


= <J = 


A = 


0, 


77 = 


-1| 


< vr TT^ 


"1 < Tn\T\OT > 


5^(00) 


= <J = 


A = 


0, 


7] = 


-1| 


< vr vr 1 


< r„|r|Or > 


Sr{ + + ) 


= <J = 


A = 


0, 


7] = 


-1| 




"1 < Tn\T\OT > 
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Generalized Bose-Einstein symmetry requires 

Sr{-+) = -^< JH <S\< T„|r|Or > (3.2) 

Using (2.4) and (2.8) we can express the amplitudes Sr{c), c = ( — h), (00), (++) in terms of 
amplitudes S!^^^ with definite total isospin / and Is 

SA-+) = -^{S^' + l^S^'} (3-3a) 

5^(00) = -^{5°° - V2520} (3.3b) 
v3 

Sr{++) = Sf (3.3c) 

Assuming that the S-matrix is invariant under the rotations in isospin space, the isospin amplitudes 
Sl^^ then do not depend on the component 1^ and we have 

Sf = Sf = Sr{++) (3.4) 

Prom (3.3a) and (3.3b) we then get the equation (1.9) 

5,(00) = S,{-+) + sl\Sr{++) (3.5) 
It is useful to write the following combinations of this equation 

f 5,(++) = 5,(00) - Sr{-+) (3.6a) 

Sr{-+) = -Sr(OO) - y|-Sr(++) (3.6b) 

5,(00) = Sr{-+) + ^Sr{++) (3.6c) 
and calculate the 5-wave intensities 

Is{c) = \Su{c)\^ + \Sd{c)\^ (3.7) 
for c = (++), ( — h), (00) using expressions on r.h.s. of (3.6). With 

Is{2) = \ls{++) (3.8) 

we then obtain 

Is{-+) + Ism - Is{'2) = +2 ^ Re[Sr{-+)S;m] (3.9a) 



Is{-+)-Ism-Is{2) = -2y|^i?e[5,(00)5*(++)] 

r 

Is{-+) - Ism + Is{2) = -2y|^i?e[5,(-+)5*(++)] 



(3.9b) 
(3.9c) 
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The interference terms on the r.h.s. of (3.9) are scalar products of four-vectors 

As{c) = {ReSu{c),ImSu{c),ReSdic), ImSd{c)} (3.10) 

in a Euchdian 4-dimensional space with the norm As{c) * As{c) = /5(c) and scalar product 

As{c)* As{c') = y%{c)./i;{d) cos n,,,{S) (3.11) 

where ^lcc'{S) is an angle between the vectors As{c) and As{c'). The relations (3.9) then read 

Is{-+) + Ism - Isi2) = +2^ Isi-+ )V Ism COS nijS) (3-12a) 
Is{-+) - Ism - Is{2) = -2v /J5(00)^/ 75(2)cosO2(5) (3.12b) 
Is{-+) - Ism + Is{2) = -2Vls{-+)y%{i)cosn3{S) (3.12c) 

The equations (3.12) represent three linearly independent constraints on the measured spectra 
Is{ — h), ^5(00) and /5(++) = §-^5(2) imposed by the requirement that the cosines have physical 
values. While the cosines are linearly independent, they satisfy a non-linear constraint 

cos^ ni{S) + cos^ n2{S) + cos^ Q'siS) - 2cosr?i(S)cosr?2(5)cosJl3(5) = 1 (3.13) 

The constraint (3.13) implies that for physical values of the cosines the phases satisfy a phase 
condition 

ni{s) + n2{s) + n^is) = (3.i4) 

There are 5 D-wave amplitudes in each process c = ( — h), (00), (++). The dominant D^{c) 
unnatural exchange amplitude corresponds to helicity A = 0. Then there are 2 unnatural 
exchange amplitudes -D^(c) and D"^ {c) and 2 natural exchange amplitudes D!^ (c) and D^^{c) 
corresponding to different combinations of helicities A = ±1 and ±2, respectively. Relations 
arising from the assumption of Generalized Bose-Einstein symmetry among L'-wave amplitudes 
for processes c = ( — h), (00), (++) are identical to relations (3.5) among the 5- wave amplitudes. 
They give rise to constraints for D-wme intensities I^jo, I^u, Idn identical in form to the 
contstraints (3.12) for 5-wave intensities. Measurements of pion pair production processes indicate 
that the helicity A = ±2 amplitudes D"^^ and D^^ are small and can be neglected at small 
momentum transfers t. While constraints analogous to (3.12) still hold, they will not be considered. 



IV. EVIDENCE FOR VIOLATION OF GENERALIZED BOSE-EINSTEIN SYMMETRY. 
A. Data used in the tests of Generalized Bose-Einstein symmetry. 

To test the constraints (3.12) we need data on the intensities Ia{c) from 7r~7r^, vr'^Tr'^ and 
ir~^7r~^ production. The available data allow to perform the tests at small momentum transfers 
0.005 < \t\ < 0.20 (GeV/c)2. 

For the 7r~7r^ channel we used two analyses of the CERN measurements on polarized target at 
17.2 GeV/c. One analysis is our high resolution analysis using Monte Carlo method presented in 
Ref. [13]. It covers mass range 580 - 1080 MeV and involves no D-wave. Another analysis is the 
CERN-Cracow-Munich (CCM) analysis ^] which covers a larger mass region of 580-1620 MeV 
and includes D-wave above 980 MeV. Both analyses are similar below 1080 MeV and produce two 
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0.005 < |t| < 0.20 (GeV/c)^ 




FIG. 1: Solutions (1,1) and (2,2) for S'-wave intensity Is{—+) for n^p vr^vr+n from high resolution 
analysis [l^ compared with Solutions 1 and 2 for intensity /s'(OO) for ir^p Tr^Tr^n from Ref. The 
intensities Is (00) were interpolated to 20 MeV bins and scaled to 17.2 GeV/c. 

solutions for the moduli \Su\'^ and \Sd\'^- Two solutions for the S'-wave intensity Is{ — h) were used 
in this mass range corresponding to combinations (1,1) and (2,2) of solutions for the moduli 

Is{-+) Solution (1,1) = |5„(1)|2 + \Sd{l)\^ (4.1) 
/5(-+) Solution (2,2) = |S42)|2 + |Srf(2)|2 

There is only one solution for the intensities Ia{ — h), A = S, , , above 980 MeV in 
the CCM analysis. The results for Is{ — h) from our analysis and the CCM analysis are shown 
in Figures 1 and 2, respectively. The results for the Z?'^-wave intensity from the CCM analysis 
are shown also in Figures 2. Figure 3 shows their results for D'-'- and D^-wave intensities. The 
intensities in the CCM analysis are in units fj,b/20 MeV. They were converted to units from 
Ref. [25}] used in our analysis using a conversion factor 0.109^6/20 MeV = 1000 events/20 MeV. 

For the vr'^vr'' channel we used the BNL data at 18.3 GeV/c [l^]. The BNL data were converted 
from native BNL units "intensity/40 MeV" into our units " 1000 events/20 MeV" using a conversion 
factor F = 0.6700 x 10~^. We obtained this factor by comparing the /2(1270) peak value in their 
Figure 5F given in units "intensity/40 MeV" with the value of coresponding 4 bins at /2(1270) 
peak in their Figure 4a given in units "events/10 MeV". The data in two bins 0.01 < \t\ < 0.10 
(GeV/c)2 and 0.10 < \t\ < 0.20 (GeV/c)^ were combined by addition to a sigle bin 0.01 < |t| < 0.20 
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0.005 < |t| < 0.20 (GeV/c)^ 




FIG. 2: 5- wave and D^-wave intensities Is{ — h) and Ido{ — ^) for tt p — > tt 7r+n from CERN-Cracow- 
Munich analysis (irt compared to two Solutions 1 and 2 for intensities /s(00) and Ido{00) for n^p 7r°7r°n 



from Ref. [2^. The intensities /s(00) and l]jo{00) were interpolated to 20 MeV bins and scaled to 17.2 
GeV/c. 

(GeV/c)^ corresponding to the CERN measurements. The data were then interpolated to 20 MeV 
bins and scaled to 17.2 GeV/c using phase and flux factor K{s,m?) given by 11] 



G{s, m 



2\ 



1 q{m?- 



Flux(s) = mp^iab 

where qijin?') = ^-s/w? — 4//^ is the pion momentum in the center of mass of the dipion system of 
mass m, and M and ^ are the nucleon and pion mass, respectively. The two Solutions 1 and 2 
for intensities /a (00), A = S,D^,D^,D^ are shown in Figures 1, 2 and 3 and compared with the 
corresponding intensities in 7r~p — > -^""vr+n. 

For the vr^vr'^ channel we used CERN data on 7r~^p — > TT~^TT~^n at 12.5 GeV/c [26i] • This data 
was used to determine the isospin 1=2 S- and D-wave amplitudes /j=o ^^'^ fj=2 ™ '""^ scattering 
using a pion exchange formula for S- and D^-wave helicity flip amplitudes 

^i(++) = Fis,t,m^)V2jTTfj=^ (4.3) 
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0.005 < |t| < 0.20 (GeV/c)^ 




FIG. 3: D^-wave and D'^-wave intensities I]ju ( — h) and I]jn {—+) for tt p 
Munich anah 



n from CERN-Cracow- 



27| compared to two Solutions 1 and 2 for intensities Ijju (00) and Ijjn (00) for tt p 



from Ref. 
GeV/c. 



The intensities liju{00) and /£)iv(00) were interpolated to 20 MeV bins and scaled to 17.2 



with an assumption that the non-flip amplitudes vanish Aq{++)=0. The partial wave intensities 
Is{++) and Ido{++) at 17.2 GeV/c were then reconstructed using 

Ia{++) = |Ai(++)|2 and Ia{2) = (4.4) 

The common factor \F{s, t, ?Ti^)p was taken from the analysis of Kaminski et al by identifying 
the mean values of our /s'(2) with the values of their 5-wave 1 = 2 contribution 1(2) at 17.2 
GeV/c and 0.005 < |t| < 0.20 (GeV/c)^ presented in their Figure 2. The errors on Ia{2) are given 
by the errors on /j=o 2 taken from the CERN analysis in Ref. [i^. The results for Is{2) 

and Ido{2) were scaled to 17.2 GeV/c and are shown in Figure 4. 

The intensities and /£)Jv(++) are not known but the CERN measurements of tt~^p 

TT^ir^n indicate they are small and consistent with zero. From the Figures 2 and 3 we see that 
Idu{ — h) and /£)]v( — h) are about 10 % of the intensity Ir)o{ — h) at the peak value for both 7r~7r+ 
and vr'^vr'^ processes. We therefore expect a similar ratio of Idu{++) and I£)n{++) to /z)o(++) 
vr+vr^" production and take 

/^^(2)=/^iv(2) = /^o(2)/10 (4.5) 
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-t^-i^ < -t < 0.20 (GeV/c)^ 




O I -I -I- - - - "I" * I I I I I I I I 

0.5 0.6 O.y 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 

mCTt^TT^) (GeV) 



FIG. 4: Isospin 1 = 2 intensities 15(2) and /£,o(2) from CERN data on 7r+7r+ production at 12.5 GeV/c [111 
scaled to 17.2 GeV/c. 

The errors on Ijju(2) and Ij^n{2) were set at 20 % of errors on l£,o{2). The sensitivity of results 
on these assumptions was tested by varying the ratios of l£)u(2) and I£)n{2) to Id°{'^) from 10 % 
to 100 %. 

B. Signatures of breakdown of Generalized Bose-Einstein symmetry. 

Generalized Bose-Einstein symmetry requires that all three cosines cosVLi.{A), k =1,2,3 for all 
amplitudes A = S, D^, and have physical values for all dipion masses for any momentum 
transfer t within the whole error volume of the measured intensities. A breakdown of the symme- 
try occurs when cosf]fc(^) calculated from constraints (3.12) and their D-waves analogues have 
unphysical values for some mass range mA_ for amplitude ^ in a large portion of the measured 
error volume. When that happens, the amplitudes violating Generalized Bose-Eistein symmetry 
no longer vanish. It is important to note that physical values of the cosines are only consistent 
with the symmetry but do not prove it. Only the unphysical values of cos Q^i^) signal conclusively 
the violation of the symmetry. However, due to possible uncertainties in the data used in the tests 
some of our results should be viewed with caution. 
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C. Test of constraints from Generalized Bose-Einstein symmetry. 

The constraints (3.12) arising from the assumption of Generalized Bose-Einstein symmetry 
were tested as follows. The error volume of the intensities Ia{ — h), -/^a(OO) and /a(2) was sampled 
by Monte Carlo method in each mass bin. For each of the 10,000 samplings used in the analysis 
equations (3.12) were used to calculate the cosines cosr2i(A), cosr22(^) and cosr23(^) for all 
solution combinations for each of the amplitudes A = S,D^,D^,D^ . A distribution of values 
of cosrifc(^) has been obtained for each k and A which defined the range and average value of 
cosr2fc(A) in each mass bin. In general, these average values of cosVtk{A)"''" were close to the mean 
values of cosVtk{A)* calculated from the mean values of the intensities. 

In each mass bin a number count was taken of the physical and unphysical values of cos^k{A) 
to quantify any possible violation of the constraints (3.12). The program also verified that the 
non-linear condition (3.13) on the cosines is satisfied for both the physical and unphysical values of 
cosVLk{A), /c=l,2,3 for each Monte Carlo sampling. Importantly, the number counts for unphysical 
value of cos Vlk{A) were identical for all three cosines in each mass bin. In the Figures 5-10 below 
we thus present only the results for cosr2i(A). 
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FIG. 5: Violations of constraints (3.12) in S'-wave pion pair production below 1080 MeV with Solution 
(1,1) for S'-wave intensity Is{ — h) from analysis [l3|. Results from CERN-Cracow-Munich analysis 271 
similar. 



The results for S-wave below 1080 MeV are shown in Figures 5 and 6. The results for S-, 
D^-, D^- and D^-waves in the mass range 980-1620 MeV are shown in Figures 7, 8, 9 and 10, 
respectively. Each Figure shows cos$7i(A) and the corresponding fraction of unphysical values of 
cosr2i(A) to evaluate the degree of violation of constraints (3.12) or, equivalently, the degree of 
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Is<00) SolTjit±on 1 



I^(OO) Solution 2 




FIG. 6: Violations of constraints (3.12) in S'-wave pion pair production below 1080 MeV with Solution 
(2,2) for S'-wave intensity Is{ — h) from analysis [l3|. Results from CERN-Cracow-Munich analysis 27 1 are 
similar. 



breaking of Generalized Bose-Einstein symmetry. 



The results for the two solutions of S'-wave intensity Is{ — h) below 1080 MeV from our high 
resolution analysis [lil] and the COM analysis 27] are very similar. Figures 5 and 6 show the 
results for the Solutions (1,1) and (2,2) of Is{ — h) from our analysis [1^]. For Solution 1 of Is(OO) 
and for both solutions of Is{ — h) there is a clear violation of the constraints (3.12) for 580 - 680 
MeV. For Solution 2 of /^(OO) there is a massive violation of constraints (3.12) in the mass range 
600-980 MeV for both solutions of /s(-+). 



Figure 7 shows the results for S'-wave in the mass range 980-1620 MeV. The constraints (3.12) 
are clearly violated at larger masses 1400-1620 MeV in both solutions for /^•(OO). Figure 8 presents 
the results for L'^-wave. Interestingly, cosQi{D^) is near 1 but clear violations of constraints 
(3.12) are again observed in the mass range 1400-1620 MeV for both solutiond of /£)o(00). 

Figures 9 and 10 present the results for D^- and D^-waves, respectively, assuming 
Ijyu{2) = Ijjn{2) is 10 % of /£)o(2). There is a large violation of consraints (3.12) in all mass 
bins for both solutions of /^(^(OO) and Idn{00). The sensitivity of this result to the assumptions 
about Iuu{2) and Idn{2) was tested by varying their ratio to lDoi2) from 10 % to 100 %. The 
large violations of constraints (3.12) persist all the way to the 100 % ratio although the fractions 
of unphysical values of cosQi{D^) and cosni{D^) get smaller at higher ratios and, at 100 %, 
become zero in some of the mass bins. 



In Figure 11 we show Solutions (1,1) and (2,2) for Is{-\ — ) from high resolution amplitude 
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FIG. 7: Violations of constraints (3.12) in S'-wave pion pair production in the mass range 980 - 1610 MeV 
with S'-wave intensity Is{ — h) from CERN-Cracow-Munich analysis 27 1. 
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FIG. 8: Violations of constraints (3.12) in ZJ^-wave pion pair production in the mass range 980 - 1610 MeV 
with Z?°-wave intensity I]jo from CERN-Cracow-Munich analysis 27 1. 
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FIG. 9: Violations of constraints (3.12) in D^-wa.ve pion pair production in the mass range 980 - 1610 MeV 
with D'^-wave intensity Ijju from CERN-Cracow-Munich analysis 27 1. 
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FIG. 10: Violations of constraints (3.12) in ZJ^-wave pion pair production in the mass range 980 - 1610 
MeV with D^-wave intensity Ijjn from CERN-Cracow-Munich analysis [27]. 
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FIG. 11: Solutions (1,1) and (2,2) for S'-wave intensity /^(H — ) in 7r+n tt+tt p at 11.85 GeV/c compared 
with Solutions 1 and 2 for 5-wave intensity /^(OO) in Tr^p ^ ir'^ir'^n at 18.3 GeV/c for momentum transfers 
0.20 < \t\ < 0.40 (GeV/c)2. 



analysis of CERN measurements of vr+n — > ■k^tt~p on polarized target at 11.85 GeV/c for dipion 
masses in the range 360 - 1040 MeV at larger momentum transfers 0.20 < |t| < 0.40 (GeV/c)^ [l^ 
and compare them with Solutions 1 and 2 for /^(OO) in the same t bin. The intensities Is{-\ — ) 
decrease rapidly with decreasing mass below 600 MeV while the Solution 1 for the intensity Is(OO) 
is large. We may assume that the intensities !${ — h) for momentum transfers 0.005 < \t\ < 0.20 
(GeV/c)^ will also rapidly decrease below 580 MeV well below the Solution 1 for the intensity 
ls(00). From this analogy we may conclude that the violation of the constraints (3.12) will remain 
and in fact increase below 580 MeV for Solution 1 of 15(00). 



D. Interpretation of the tests. 

Despite some caveats regarding the accuracy of the data used in the tests, the tests reveal a 
rather convincing case for violation of Generalized Bose-Einstein symmetry and for the presence 
of amplitudes violating this symmetry in 'ir~p — > ■K~'K~^n reaction. The violation of the generalized 
Bose-Einstein symmetry is not related to strong interactions or to electromagnetic interactions. 
It arises from a new CPT violating interaction of the pion creation process with a quantum 
environment. The presence of p'^(770) in the 5- wave due to this interaction is associated with 
violations of the constraints by the Solution 1 for 15(00) below 680 MeV and by Solution 2 for 
15(00) for masses below 980 MeV. 

The mass range of 1400 - 1620 MeV is within the mass region of ^"(1450) resonance. The 
violations of constraints (3.12) in this mass range by 5- and D'^-wave amplitudes suggests the 
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presence of p°(1450) in the 5- and D^-waves akin to p^{770) — /o(980) mixing. This result is 
expected from our model of interaction of the pion creation process with the environment [l^ 
in which the interaction with the environment induces transitions of p^{lA50) into the S- and 
-waves. There may be a similar effect in the D^- and D^- waves. 



V. EVIDENCE FOR QUANTUM ENTANGLEMENT OF tt 7r+ ISOSPIN STATES. 

When Generalized Bose-Einstein symmetry is violated the two-pion state tt^tt^ produced in 
TT~p — > Tr~7T^n needs no longer be a separable state (2.3). In general, it will be an entangled state 

|£'(7r"7r+) >= as\S > +aA\A >= a|7r"7r+ > +6|7r+7r" > (5.1) 

where IS* > and \A > are symmetric and antisymmetric tt^tt^ isospin states (2.4), respectively, 
and 

a = -j^ias + aA) (5.2) 
b = -^{as - Qa) 

Only for ac = = A= are the two-pion states 7r"7r+ separable. The transversity amplitudes then 

v2 

have a general form 

J=even: A^"^^ = asCi^ + gaV^^^ (5.3) 

J=odd: Ai"^^^ = aACil + asV^^l 

where C^'J. and are Generalized Bose-Einstein symmetry conserving and violating amplitudes. 
The entanglement amplitudes as and a a do not depend on transversity t = u,d nor on solution 
labels i,j = 1,2. We can see that from the general expression for helicity amplitudes in terms of 
transversity amplitudes [l^, [3] 

= + (-^r<'M (5-4) 

where n = 0, 1 for helicity non-flip and flip amplitudes, respectively. The helicity amplitudes 
A'^^{ij) must share the same isospin state (5.1) as the transversity amplitudes in order for them 

to be expressed in terms of symmetry conserving and violating helicity amplitudes C^f^ and 

formed by combinations (5.4) of symmetry conserving and violating transversity amplitudes C^'J. 

and V^/^, respectively. 



To examine the data below 1080 MeV for evidence of quantum entanglement of vr vr^ isospin 
states we shall focus on S- and P-wave amplitudes with helicity A = which we write in the form 

Sr{i) = asS^{i) + aAS]:{i) (5.5) 
Lrii) = aAL^{i) + asL^ii) 

where i = 1,2. The symmetry conserving amplitude S^{i) must satisfy relation (3.5) 

5f (i) = Sr{00, i) - y|^r(++, i) (5.6) 
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where Sr{00,i) and Sr{++,i), i=l,2 are transversity S'-wave amplitudes in 7r~p — > vr'^vr'^n and 
vr+p 7r~^7r~^n reactions, respectively. From our model for non- unitary dynamics of p^{770) — 
/o(980) mixing [l^ we have for the Generalized Bose-Einstein symmetry violationg amplitudes 

S^{i) ~ L^{^) (5.7) 
lU^) - 5f (^) 

where S!^^{i) is the isospin 7 = component of S'^(i). We are interested in the difference 

5.(1) - 5.(2) = asiS^il) - 5f (2)) + a^(5]f (1) - 5^(2)) (5.8) 

To proceed we shall make several plausible assumptions. First, we expect that there are very small 
differences between the 1 = 2 amplitudes (which themselves are small below ~ 800 MeV) and set 

5.(++,l) -5.(++,2) = (5.9) 

Then 

5f (1) - 5f (2) = 5.(00, 1) - 5.(00, 2) (5.10) 

From the data in Figure 1 we expect this difference to be large below ~ 900 MeV. In our next 
step we note that the amplitudes |L.(l)p and |L.(2)p are very similar for both transversities 
T = u,d [13]. We thus expect there will be only small differences between the amplitudes 5)^(1) 
and 5)^(2) and set 

5^(1) - 5^(2) = (5.11) 

From (5.8) we thus have 

|5.(1) - 5.(2)|2 = |asp|5.(00,l) - 5.(00, 2)1^ (5.12) 

In our previous work [3] we have shown that for physical helicity amplitudes the phases of the 
two solutions 5.(1) and 5.(2) must be equal. Then 

|5.(1) - 5.(2)|2 = (|5.(1)| - |5.(2)|)2 (5.13) 

This term is known and the moduli squared of the amplitudes are shown in Figure 12 taken 
from Since we do not know the relative phase of the amplitudes 5.(00,1) and 5.(00,2) we 
rewrite (5.12) as inequalities 

|asP(|5.(00,l)| - 15.(00,2)1)2 < (|^^(^)| _ \s^^2)\f (5.14) 
(|5.(1)| -|5.(2)|)2 < |asP(|5.(00, 1)1 + 15.(00,2)1)2 

In our final step we assume that 

|5„(00, 1)|2 = |5rf(00, 1)|2 = /s(00, l)/2 (5.15) 
|5„(00,2)|2 = |5rf(00,2)|2 = /s(00,2)/2 

where /s(00, 1) and /5(00, 2) are the Solutions 1 and 2 for the 5-wave intensity in Tr~p — > Tr^n^n 
shown in Figure 1. This assumption is in fact an assumption made in BNL data analysis 22l |. 
From (5.14) we finally obtain lower and upper bounds on the entanglement amplitude \as\ 



V2 ll^-WI- 1^.(2)11 < i^^i < ^ 11^.(1)1-1^.(2 



V/s(00,l) + V/5(00,2) V^5(00,l)- V^5(00,2) 



(5.16) 
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FIG. 12: Solutions 1 and 2 for 5-wave amplitudes l^up and jS'dp in tt p ^ tt 7r+n at 17.2 GeV/c from 
analysis [l^ . 
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FIG. 13: Estimates of lower and upper bounds on entanglement amplitudes \as\ and |a^| from data on 
n^p TT^TT+n and ir^p — > Tr^Tr'^n at 17.2 GeV/c. 
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From the normalization condition \as\'^ + joyip = 1 we obtain lower and upper bounds on the 
entanglement amplitude IqaI- 

The results for entanglement amplitudes \as\ and {uaI calculated for both cases of transversity 
T = u,d are shown in Figure 13. To estimate \as\ we used mean values of |5'T-(i)| and 15(00, i), 
i = 1,2 and errors were not calculated. Below 840 MeV the results from both transversity 
calculations are consisten with \as\ ~ 0.10 — 0.20 and |a^| ~ 0.98 — 0.99 indicating a quantum 
entanglement of the vr^vr^ isospin states in tt~p — > 7r~7r~*'n. 

The effect of the quantum entanglement below 840 MeV is a suppression of the conserving 
amplitude S^{i) and the clear dominance of the violating amplitude {i) in the mass range of 
p°(770) resonance in the S'-wave amplitudes Sr{i)- For dipion masses 980-1080 MeV calculations 
with both transversities are again consistent and suggest larger values of \as\- This value of 
\as\ leads to large contribution from the violating amplitude L^ii) to the amplitudes Lt-(0 
which manifests itself by the presence of a dip in at ~980 MeV corresponding to /o(980) 

resonance. There is a divergence between the two calculations in the mass range 840-980 MeV 
which suggests that some of our simplifying assumptions are violated. This violation may be also 
responsible for the large range of values of the entanglement amplitudes above 980 MeV. 



We cannot estimate the entanglement amplitudes above 980 MeV using the CERN-Cracow- 



Munich analysis 27|. The S- and P-wave subsystem of the reduced density matrix measured 



on transversely polarized target is always analytically solvable and yields two solutions for each 



transversity amplitude Au{i), Ad{j),i, j = 1,2 [13|. While CERN-Cracow-Munich analysis used 



the analytical solutions below 980 MeV, above this mass fits were made to measured moments 
which are linear combinations of density matrix elements. These fits produced a single solution for 
the amplitudes. This can happen if the analytical solutions are close and cannot be distinguished 
by the fits. The single solution for the S'-wave amplitudes obtained in the fits cannot be used in 
(5.16) to estimate the entanglement amplitudes. 



VI. CONCLUSIONS. 



We have derived three linearly independent constraints on partial wave intensities with even 
spin in TT~p — > 7r~7r"'"n, 7r~p — > vr'^vr'^n and vr+p — > vr^vr+n reactions that must be satisfied by the 
data if the Generalized Bose-Einstein symmetry holds true. Available data violate the constraints. 
The violation of the symmetry implies the presence of the symmetry violating contributions to 
transversity amplitudes in 'ir~p — > 7r~7r"'"n and allows for quantum entanglement of the tt^tt^ 
isospin states. We have derived approximate expressions for lower and upper bounds on entan- 
glement amplitudes. The bounds provide clear experimental evidence for quantum entanglement 
of vr^TT^ isospin states below 840 MeV and suggest the entanglement at higher dipion masses. 
The small values of \as\ ~ 0.10 — 0.20 below 840 MeV explain the puzzling difference between 
the S-wave intensities in ir^p ■K~TT~^n and 7r~p — > ir^ir^n and reveal a suppression of isospin 
/ = 0, 2 contribution in the S'-wave amplitudes in ir^p vr^vr"^?!. The large isospin 1=1 
contribution of p^(770) to both S- and P-wave amplitudes is due to large entanglement amplitude 
\aA\ ~ 0.98 - 0.99. 

A crucial aspect of the our tests was the treatment of all solutions for the intensities as valid 
physical solutions. For instance, when looking at the tests of the constraints (3.12) for the 5- wave 
amplitudes below 1080 MeV it is not enough to look only at the Solution 1 for /s'(OO) or to look at 
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the two solutions for Is (00) separately. Instead, the tests for the two solutions must be viewed and 
interpreted together since the underlying amplitudes for different solutions form a single system. 

The violations of Generalized Bose-Einstein symmetry and the entanglement of 7r^7r+ isospin 
states are not related to strong or electromagnetic interactions. We suggest that they provide new 
evidence for the existence of a quantum environment and its CPT violating interaction with pion 
creation processes. The origin of the environment remains an open question. 



The new interaction manifests itself in several ways. It leads to splitting of the transversity 
amplitudes and to two physical solutions Au{i),i = 1,2 and Aii{j),j = 1,2 which 
results in mixed final state (1.4) and in a CPT violating non-unitary evolution of pure initial 
states to mixed final states in all vrA^ ttttN processes 15, The splitting reflects the inter- 
acting degrees, or the quantum numbers, of the environment described by quantum states \i > \j >. 



The interaction also leads to isospin conserving transitions between resonant qq modes with 
spin K, helicity /i and isospin Ik produced in n^p — > qqn and two-pion states with spin J, 
helicity A and isospin 1(11^) = Ik- As the result of these transitions resonant modes can "leak" 
into different partial wave amplitudes The transitions are observable as violations of 

Generalized Bose-Einstein symmetry, p^(770) — /o(980) mixing [131 and the presence of ^"^(1450) 
in S- and i^'^-wave amplitudes. Transitions from K = 2 qq mode to J = vr'^vr'^ state could also 
explain the large differences in the two solutions for the S'-wave amplitudes in 7r~p — > vr'^vr'^n below 
900 MeV. 



The CPT violating interaction with the environment produces entangled two-pion isospin state 
|£'(7r^7r^) >= as\S > +aA\A >. A separable state tt^tt^ is expected from strong interactions 
and Generalized Bose-Einstein symmetry. The entanglement of the 7r~7r"'" state is another aspect 
of CPT violation in tt~p — > 7r~7r+n. The entangled state is not experimentally preparable since 
the entanglement amplitudes are not constants but may depend on energy s, dipion mass m 
and momentum transfer t. As the result the amplitude for the CPT conjugate process is not 
defined and CPT symmetry looses its meaning in the interaction of pion creation process with 
the quantum environment. 

Our study has been limitted to partial wave intensities since the reactions iT~p vr^vr'^n and 
TT~^p 7r+7r"^n were measured only on unpolarized targets. Future experiments may measure all 
pion creation processes on polarized targets and thus enable more detailed studies of quantum 
entanglement and non- unitary dynamics of pion creation processes on the level of amplitudes, 
deepening our understanding of this new class of phenomena. 
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